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Description 



[DIRECTIONAL ELECTROMAGNETIC 
WAVE RESISTIVITY APPARATUS AND 
METHOD] 

Cross Reference to Related Applications 

[0001] This application claims the priority from U.S. Provisional 
Application No. 60/472,686 filed on May 22, 2003. This 
Provisional Application is incorporated by reference in its 
entirety. 
Background of Invention 

[0002] Field of the Invention 

[0003] The invention relates generally to the field of well logging. 
More particularly, the invention relates to improved tech- 
niques in which instruments equipped with antenna sys- 
tems having transverse or tilted magnetic dipole repre- 
sentations are used for electromagnetic measurements of 
subsurface formations and for placing wells with respect 
to geological boundaries in a reservoir. The invention has 



general application in the well logging art, but is particu- 
larly useful in logging-while-drilling operation. 
[0004] Background of the Related Art 

[0005] Various well logging techniques are known in the field of 
hydrocarbon exploration and production. These tech- 
niques typically use instruments or tools equipped with 
sources adapted to emit energy into a subsurface forma- 
tion that has been penetrated by a borehole. In this de- 
scription, "instrument" and "tool" will be used inter- 
changeably to indicate, for example, an electromagnetic 
instrument (or tool), a wire-line tool (or instrument), or a 
logging-while-drilling tool (or instrument). The emitted 
energy interacts with the surrounding formation to pro- 
duce signals that are then detected and measured by one 
or more sensors. By processing the detected signal data, a 
profile of the formation properties is obtained. 

[0006] Electromagnetic (EM) induction and propagation logging 
are well-known techniques. The logging instruments are 
disposed within a borehole to measure the electrical con- 
ductivity (or its inverse, resistivity) of earth formations 
surrounding the borehole. In the present description, any 
reference to conductivity is intended to encompass its in- 
verse, resistivity, or vice versa. A typical electromagnetic 



resistivity tool comprises a transmitter antenna and one or 
more (typically a pair) receiver antennas disposed at a dis- 
tance from the transmitter antenna along the axis of the 
tool (see Figure 1). 
[0007] induction tools measure the resistivity (or conductivity) of 
the formation by measuring the voltage induced in the re- 
ceiver antenna(s) as a result of magnetic flux induced by 
currents flowing through the emitting (or transmitter) an- 
tenna. An EM propagation tool operates in a similar fash- 
ion but typically at higher frequencies than do induction 
tools for comparable antenna spacings (about 10 6 Hz for 
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propagation tools as compared with about 10 Hz for the 
induction tools). A typical propagation tool may operate at 
a frequency range of 1 kHz - 2 MHz. 
[0008] Conventional transmitters and receivers are antennas 

formed from coils comprised of one or more turns of in- 
sulated conductor wire wound around a support. These 
antennas are typically operable as sources and/or re- 
ceivers. Those skilled in the art will appreciate that the 
same antenna may be use as a transmitter at one time and 
as a receiver at another. It will also be appreciated that the 
transmitter-receiver configurations disclosed herein are 
interchangeable due to the principle of reciprocity, i.e., 



the "transmitter" may be used as a "receiver", and vice- 
versa. 

[0009] The antennas operate on the principle that a coil carrying 
a current (e.g., a transmitter coil) generates a magnetic 
field. The electromagnetic energy from the transmitter 
antenna is transmitted into the surrounding formation, 
and this transmission induces eddy currents flowing in the 
formation around the transmitter (see Figure 2A). The 
eddy currents induced in the formation, which are func- 
tions of the formation's resistivity, generate a magnetic 
field that in turn induces an electrical voltage in the re- 
ceiver antennas. If a pair of spaced-apart receivers is 
used, the induced voltages in the two receiver antennas 
will have different phases and amplitudes due to geomet- 
ric spreading and absorption by the surrounding forma- 
tion. The phase difference (phase shift, O and amplitude 
ratio (attenuation, A) from the two receivers can be used 
to derive the resistivity of the formation. The detected 
phase shift (O) and attenuation (A) depend on not only the 
spacing between the two receivers and the distances be- 
tween the transmitter and the receivers, but also the fre- 
quency of EM waves generated by the transmitter. 

[0010] | n conventional induction and propagation logging instru- 



merits, the transmitter and receiver antennas are mounted 
with their axes along the longitudinal axis of the instru- 
ment. Thus, these tools are implemented with antennas 
having longitudinal magnetic dipole (LMD) representa- 
tions. An emerging technique in the field of well logging 
is the use of instruments including antennas having tilted 
or transverse coils, i.e., where the coil's axis is not parallel 
to the longitudinal axis of the tool. These instruments are 
thus implemented with a transverse or tilted magnetic 
dipole (TMD) antenna. Those skilled in the art will appre- 
ciate that various ways are available to tilt or skew an an- 
tenna. Logging instruments equipped with TMD antennas 
are described, e.g., in: U.S. Patent Nos. 6,163,155; 
6,147,496; 5,115,198; 4,319,191; 5,508,616; 5,757,191; 
5,781,436; 6,044,325; and 6,147,496. 
1 ] Figure 2A presents a simplified representation of eddy 
currents and electromagnetic (EM) energy flowing from a 
logging instrument disposed in a borehole portion or seg- 
ment that penetrates a subsurface formation in a direction 
perpendicular to the sedimentation layers. This is not, 
however, an accurate depiction of all the numerous seg- 
ments that make up a borehole particularly when the 
borehole has been directionally-drilled as described be- 



low. Thus, segments of a borehole often penetrate forma- 
tion layers at an angle other than 90 degrees, as shown in 
Figure 2B. When this happens, the formation plane is said 
to have a relative dip. A relative dip angle, O, is defined as 
the angle between the borehole axis (tool axis) BA and the 
normal N to the plane P of a formation layer of interest. 

[0012] Drilling techniques known in the art include drilling bore- 
holes from a selected geographic position at the earth's 
surface, along a selected trajectory. The trajectory may 
extend to other selected geographic positions at particu- 
lar depths within the borehole. These techniques are 
known collectively as "directional drilling" techniques. One 
application of directional drilling is the drilling of highly 
deviated (with respect to vertical), or even horizontal, 
boreholes within and along relatively thin hydrocarbon-bear- 
ing earth formations (called "pay zones") over extended 
distances. These highly deviated boreholes are intended 
to greatly increase the hydrocarbon drainage from the pay 
zone as compared to "conventional" boreholes which "ver- 
tically" (substantially perpendicularly to the layering of the 
formation, as shown in Figure 2A) penetrate the pay zone. 

[0013] | n highly deviated or horizontal borehole drilling within a 
pay zone, it is important to maintain the trajectory of the 



borehole so that it remains within a particular position in 
the pay zone. Directional drilling systems are well known 
in the art which use "mud motors" and "bent subs," as well 
as other means, for controlling the trajectory of a bore- 
hole with respect to geographic references, such as mag- 
netic north, the earth's gravity (vertical), and the earth's 
rotational velocity (with respect to inertial space). Layering 
of the formations, however, may be such that the pay 
zone does not lie along a predictable trajectory at geo- 
graphic positions distant from the surface location of the 
borehole. Typically the borehole operator uses informa- 
tion (such as LWD logs) obtained during borehole drilling 
to maintain the trajectory of the borehole within the pay 
zone, and to further verify that the borehole is, in fact, 
being drilled within the pay zone. 
[0014] Techniques known in the art for maintaining trajectory are 
described for example in Tribe et al., Precise Well Placement 
using Rotary Steer able Systems and LWD Measurement, Society of 
Petroleum Engineers, Paper 71396, September 30, 2001. 
The technique described in this reference is based upon 
LWD conductivity sensor responses. If, as an example, the 
conductivity of the pay zone is known prior to penetration 
by the borehole, and if the conductivities of overlying and 



underlying zones provide a significant contrast with re- 
spect to the pay zone, a measure of formation conductiv- 
ity made while drilling can be used as a criterion for 
"steering" the borehole to remain within the pay zone. 
More specifically, if the measured conductivity deviates 
significantly from the conductivity of the pay zone, this is 
an indication that the borehole is approaching, or has 
even penetrated, the interface of the overlying or underly- 
ing earth formation. As an example, the conductivity of an 
oil-saturated sand may be significantly lower than that of 
a typical overlying and underlying shale. An indication 
that the conductivity adjacent the borehole is increasing 
can be interpreted to mean that the borehole is approach- 
ing the overlying or the underlying formation layer (shale 
in this example). The technique of directional drilling us- 
ing a formation property measurement as a guide to tra- 
jectory adjustment is generally referred to as "geosteer- 
ing." 

[0015] | n addition to EM measurements, acoustic and radioactive 
measurements are also used as means for geosteering. 
Again using the example of an oil-producing zone with 
overlying and underlying shale, natural gamma radioactiv- 
ity in the pay zone is generally considerably less than the 



natural gamma ray activity of the shale formations above 
and below the pay zone. As a result, an increase in the 
measured natural gamma ray activity from a LWD gamma 
ray sensor will indicate that the borehole is deviating from 
the center of the pay zone and is approaching or even 
penetrating either the upper or lower shale interface. 

[0016] |f f as j n the prior examples, the conductivity and natural 
radioactivity of the overlying and underlying shale forma- 
tions are similar to each other, the previously described 
geosteering techniques indicate only that the borehole is 
leaving the pay zone, but do not indicate whether the 
borehole is diverting out of the pay zone through the top 
of the zone or through the bottom of the zone. This 
presents a problem to the drilling operator, who must 
correct the borehole trajectory to maintain the selected 
position in the pay zone. 

[0017] em induction logging instruments are well suited for 
geosteering applications because their lateral (radial) 
depth of investigation into the formations surrounding the 
borehole is relatively large, especially when compared to 
nuclear instruments. The deeper radial investigation en- 
ables induction instruments to "see" a significant lateral 
(or radial) distance from axis of the borehole. In geosteer- 



ing applications, this larger depth of investigation enables 
the detection of approaching formation layer boundaries 
at greater lateral distances from the borehole, which pro- 
vides the drilling operator additional time to make any 
necessary trajectory corrections. Conventional propaga- 
tion-type instruments are capable of resolving axial and 
lateral (radial) variations in conductivity of the formations 
surrounding the instrument, but the response of these in- 
struments generally cannot resolve azimuthal variations in 
the conductivity of the formations surrounding the instru- 
ment. Furthermore, such instruments are unable to sense 
anisotropy in vertical wells. 
[0018] t wo important emerging markets make the removal of 

these shortcomings more urgent. The first emerging field 
is the increasing need for accurate well placement, which 
requires directional measurements to make steering deci- 
sions to place the borehole optimally in the reservoir. The 
second is the low resistivity pay in laminated formations 
where accurate identification and characterization of hy- 
drocarbon reserves is not possible without knowing the 
resistivity anisotropy. Many recent patents disclose meth- 
ods and apparatus to make directional measurements and 
obtain resistivity anisotropy. For logging while drilling ap- 



plications, U.S. Patent No. 5,508,616 to Sato et al. dis- 
closes an induction-type tool with two coils tilted at dif- 
ferent directions not aligned with the tool's longitudinal 
axis. The directionality of the measurement is illustrated 
through a simple argument that the sensitivity function of 
the two tilted coils is concentrated towards the overlap- 
ping region of the sensitivity area of each coil. Through 
rotation of the tool, Sato et al claim that a deep azimuthal 
resistivity image of the formation can be obtained. How- 
ever, this patent reference does not provide any details as 
to how the azimuthal resistivity can be obtained, nor does 
it describe any further boundary detection/charac- 
terization techniques required for quantitative geosteering 
decision-making. 
[° 019 ] U.S. Patent No. 6,181,138 to Hagiwara and Song extends 
Sato et al's single fixed directional coils into co-located 
triple orthogonal induction coils at the transmitter and re- 
ceiver location. No tool rotation is said to be required, 
since the focusing direction can be tuned to arbitrary ori- 
entation through linear combination of the orthogonal coil 
responses. It is not clear if there is a shield design that 
will allow the passing of all the required EM components 
without severe uncontrollable distortion of the wave form 



for "while drilling" applications. 
[0020] u.S. Patent No. 6,297,639 to Clark et al., assigned to the 
assignee of the present invention, discloses methods and 
apparatus for making directional measurements utilizing 
various shield designs to provide selected attenuation of 
EM wave energy for axial, tilted, and transverse antenna 
coils. This patent reference describes, among other 
things, general directional induction and propagation 
measurements with tilted coils and appropriate shields, 
along with a process for conducting borehole compensa- 
tion for them which is non-trivial. A one-axial and one- 
tilted transmitter/receiver coil combination is explicitly 
described by Clark et al., along with its application for bed 
boundary direction detection by observing azimuthal vari- 
ation of the induced signal as the tool rotates. The az- 
imuthal variation of the coupling can be used for steering 
wells while drilling. More shield patents have since been 
granted, including U.S. Patent No. 6,351,127 to Rosthal et 
al., and U.S. Patent No. 6,566,881 to Omeragic et al, both 
of which are assigned to the assignee of the present in- 
vention. 

[0021] u.S. Patent No. 6,476,609 to Bittar extends an earlier 
anisotropy patent describing both transmitters and re- 



ceivers possibly having a tilt angle , U.S. Patent No. 
6,163,155 also to Bittar, to the area of geosteering appli- 
cation. The bedding response of up/down tilted induction 
and propagation apparatus is described through the dif- 
ference or ratio of signals at two different orientations, 
but no shielding is mentioned. Nor are the effects of 
anisotropy or dipping considered. Also lacking is a de- 
scription of how to use these measurements to derive a 
precise distance to a formation bed boundary. The "609 
patent implicitly assumes that bedding orientation is pre- 
cisely known so as to calculate the up/down response. No 
technique, however, is disclosed to locate the precise up 
or down direction prior the calculation of the up-down di- 
rectional signals. 
[0022] U.S. Patent Application Publication No. 2003/0085707 to 
Minerbo et al, assigned to the assignee of the present in- 
vention, discloses tool configurations and symmetrization 
techniques that simplify the response of the directional 
measurements to the point that it becomes almost inde- 
pendent to anisotropy or dipping angle. Responses to bed 
boundary distance with different dip and anisotropy es- 
sentially overlap except near the bed boundary. Both two- 
coil (one transmitter and one receiver: "TR") induction 



style and three-coil (one transmitter and two receivers: 
"TRR") propagation-style measurements can be sym- 
metrized to achieve this simplification. The symmetriza- 
tion is done between two tilted TR pairs of the same spac- 
ing, but with the transmitter tilted angle and receiver 
tilted angle exchanged. Only cases where the magnetic 
moments of the transmitters and receivers are lying in the 
same plane is considered. This has a disadvantage of not 
being able to provide the required signal for geosteering 
all the time during sliding, which is the case for well 
placement with a mud motor during the angle build on 
the trajectory. If the magnetic moment of the tool hap- 
pens to lie parallel to the bedding during sliding, the up/ 
down directional signal generated will be zero indepen- 
dent of the distance to the boundary. Thus no monitoring 
of the distance to the boundary is possible. 
[0023] U.S. Patent Application Publication No. 2003/0200029 to 
Omeragic et al, also assigned to the assignee of the 
present invention, discloses propagation-style directional 
measurements for anisotropy determination in near- 
vertical wells with borehole compensation. Inversion tech- 
niques are also used to obtain the anisotropic formation 
property. U.S. Patent Application Publication No. 



2003/0184302 to Omeragic and Esmersoy, assigned to 
the assignee of the present invention, also discloses tech- 
niques for looking-ahead with directional measurements. 

[00 24 ] U.S. Patent Application Publications No. 2004/0046560A1 
and 2004/0046561A1 to Itzkovicz et al, discloses the use 
of quadrupole antennas, and transverse dipole- 
quadrupole coupling and induction style measurements 
with similar directional characteristics to a conventional 
cross-dipole XZ response. Practical realization on a 
metallic collar and adequate shielding of such antennas is 
not clear. Also, the borehole effect of such measurements 
and its interaction/coupling with the boundary effect may 
be different from XZ-style measurements. 

[0025] None of the above patent references discloses the use of 
detailed azimuthal responses of the measured signal or 
techniques to extract such responses. These references 
further fail to teach how to use the directional measure- 
ment to arrive at boundary distances for geosteering use. 
Only the so-called up/down measurement, which is the 
difference in the measured signal between the tool focus- 
ing directly towards and away from the formation bed, is 
mentioned. The precise bedding dip and azimuth infor- 
mation is usually not known before the drilling, and they 



also frequently vary in challenging well placement situa- 
tions where geosteering is required. Using a predefined 
bedding up/down direction produces at best degraded 
measurement and at worst can lead to wrong geosteering 
decisions when the bedding azimuth suddenly changes. In 
principle, the measurements can be binned azimuthally 
downhole. This technique has a number of drawbacks in- 
cluding difficulties in aligning the top and bottom bins 
precisely with the orientation of the formation bedding, 
and the inability to use (i.e., wasting) the data that are not 
in the up and down bins. The large memory required to 
record the azimuthal data with sufficient accuracy is also 
an issue. 

[0026] More importantly, the existing art for geosteering using 
directional measurements works only for steering up and 
down. There are many cases where the wellbore has to 
move azimuthally to avoid exiting the pay zone. 

[0027] a need therefore exists for methods and techniques of 

extracting and analyzing the azimuthal dependence of di- 
rectional logging measurements, using measurements 
taken at all the azimuthal angles, for characterizing the 
earth formation and for steering wells during drilling with 
improved accuracy. 



[0028] a need further exists for providing the bedding azimuth 
from the directional measurements, and generating mea- 
surements that can be used for well placement in up/ 
down or azimuthal steering. 

[0029] a need further exists for methods of utilizing these direc- 
tional measurements in real-time to obtain bed boundary 
distances and to obtain accurate earth models such that a 
geosteering decisions can be made for well placement. 

[0030] a need further exists for a method of detecting the pres- 
ence of resistivity anisotropy in formation layers adjacent 
near-vertical wells. 

[0031] a still further need exists for an efficient system that pro- 
vides such directional measurements, analyzes them 
downhole, and transmits relevant information to surface 
to facilitate geosteering up/down or azimuthally during 
well placement. It would be further advantageous if such a 
system could provide distance-to-boundary information 
during the sliding phases of drilling (i.e., no drill string 
rotation) as well as when the system / tool is rotating. 

[0032] DEFINITIONS 

[0033] Certain terms are defined throughout this description as 
they are first used, while certain other terms used in this 
description are defined below: 



[0034] "Convergence" means the condition when iteratively-cal- 
culated values approach observed values or finite limits as 
the number of iteration cycles increases. 

[0035] "Crossplot" means a graph, or the creation of such a 

graph, that indicates the relationship between two differ- 
ent measurements of the same subject or sample. 

[0036] "inversion" or "invert" means deriving a model (a.k.a. "in- 
version model") from measured data (e.g., logging data) 
that describes a subsurface formation and is consistent 
with the measured data. 

[0037] "Toolface" refers to the angular orientation of an instru- 
ment about its longitudinal axis, and represents an angle 
subtended between a selected reference on the instru- 
ment's housing (e.g., a drill collar) and either the gravita- 
tionally uppermost wall of the wellbore or geographic 
north. 

[0038] "Symmetry" or "symmetric," as used herein, refers to a 
configuration in which sets of transmitter-receiver ar- 
rangements are provided in opposite orientations along 
the longitudinal axis of a tool, such that these transmit- 
ter-receiver sets can be correlated with a standard sym- 
metry operation (e.g., translation, mirror plane, inversion, 
and rotation) with respect to a point on the tool axis or a 



symmetry plane perpendicular to the tool axis. 
Summary of Invention 

[0039] | n one aspect, the present invention provides a novel on- 
the-fly data processing technique to extract signals from 
the azimuthal variation of the directional measurement 
that are relevant for formation characterization and 
geosteering use. Instead of placing logging data into az- 
imuthal bins and then determining the values associated 
with up and down directions as is done in conventional 
imaging and proposed by others, the present invention 
takes advantage of the simplicity of the physics of the 
logging response. More particularly, the relevant bound- 
ary, anisotropy and fracture signals are extracted from the 
formation response through fitting of the azimuthal varia- 
tion of the measured voltages to some sinusoidal func- 
tions. The orientation of the bedding is also obtained as a 
result. This on-the-fly processing improves the accuracy 
of the measurements because data in all azimuthal direc- 
tions, or angles, are used and no quantization according 
to particular angles occurs. Such large processing is made 
possible with integer calculations in a digital signal pro- 
cessor (DSP), which also represents an innovation in the 
implementation. 



[0040] By extracting the relevant voltage couplings according to 
their order of azimuthal dependence, the present inven- 
tion allows for calculations of propagation style measure- 
ments with only one transmitter and receive pair. It also 
allows for superposition of measurements of different TR 
pairs of even different azimuthal orientation to generate 
other measurements of unique property, such as the sym- 
metrized or anti-symmetrized measurements. 

[0041] The present invention further provides an innovative mea- 
surement sequence to increase the tolerance to irregular 
rotation, harmonic binning, and even stick-slip. A rapid 
firing sequence is utilized, and this sequence is random- 
ized at each acquisition cycle to reduce cyclic lock-up at 
particular rotational speeds. 

[0042] Accordingly, the present invention may be expressed as a 
method for characterizing a subsurface formation, begin- 
ning with the step of disposing a suitable logging instru- 
ment within a borehole. The logging instrument is 
equipped with at least first transmitter and receiver an- 
tennas spaced apart by a first distance. At least one of the 
first antennas has a tilted magnetic dipole with respect to 
the longitudinal axis of the instrument. The first antennas 
are oriented about the axis of the logging instrument such 



that the at least one tilted magnetic dipole corresponds to 
a first azimuthal angle. The logging instrument is az- 
imuthally-rotated within the borehole, e.g., by rotation of 
the drill collar or drill string tool containing the instru- 
ment. While the logging instrument is rotating, the first 
transmitter antenna is activated to transmit electromag- 
netic energy into the formation. Also, while the logging 
instrument is rotating, a set of first voltage signals associ- 
ated with the transmitted electromagnetic energy using 
the first receiver antenna is directionally measured as a 
function of the azimuthal orientation of the logging in- 
strument. The directional measurements determine the 
azimuthal variation of the measured first voltage signals. 
This azimuthal variation is fitted to approximate func- 
tions. The activating, measuring, and fitting steps may be 
repeated to execute subsequent acquisition cycles. 

[0043] | n a particular embodiment, the fitting step is executed 

while the first voltage signals are being measured and the 
fitting is stopped when convergence has been achieved. 
The fitting coefficients are preferably determined using a 
Fast Fourier Transform. 

[0044] | n a particular embodiment, the fitting functions are sinu- 
soids defined by coupling components of the first trans- 



mitter antenna's magnetic dipole and first receiver an- 
tenna's orientation vectors. The coefficients of the fitting 
components are preferably functions of earth formation 
parameters including at least one of resistivity of forma- 
tion beds, location of the logging instrument, borehole 
deviation, azimuth angle at the location of the logging in- 
strument, and a combination thereof. The fitting coeffi- 
cients preferably include constant, sinO, cosO, sin 24> and 
cos20 terms that define an iterative fitting algorithm use- 
ful for determining the azimuthal dependence of the di- 
rectional measurements. 
[0045] As mentioned above, the present invention is adaptive to 
the superposition of measurements of different transmit- 
ter-receiver ("TR") pairs. Accordingly, in a particular em- 
bodiment, the logging instrument is further equipped with 
second transmitter antenna and receiver antennas spaced 
apart by the first distance. The second transmitter has a 
magnetic dipole whose tilt corresponds to the tilt of the 
first receiver antenna, and the second receiver antenna 
has a magnetic dipole whose tilt corresponds to the tilt of 
the first transmitter antenna, such that at least one of the 
second antennas has a tilted magnetic dipole. The second 
transmitter and receiver antennas are oriented about the 



axis of the logging instrument such that the at least one 
tilted magnetic dipole corresponds to a second azimuthal 
angle. Thus, while the logging instrument is rotating, the 
second transmitter antenna to transmit electromagnetic 
energy into the formation, and a second set of voltage 
signals associated with the transmitted electromagnetic 
energy are directionally measured using the second re- 
ceiver antenna, as a function of the azimuthal orientation 
of the logging instrument. These directional measure- 
ments determine the azimuthal variation of the measured 
second voltage signals. As with the measured first voltage 
signals, the azimuthal variation of the measured second 
voltage signals are fitted to approximate functions. 

[0046] | n a particular embodiment, the second azimuthal angle 
differs from the first azimuthal angle by substantially 90 
degrees. Alternatively, the second azimuthal angle may be 
substantially equal to the first azimuthal angle. 

[0047] | n a particular embodiment, the fitting functions are sinu- 
soids defined by coupling components of the first trans- 
mitter antenna's magnetic dipoles and first receiver an- 
tenna's orientation vectors, and by coupling components 
of the second transmitter antenna's magnetic dipoles and 
second receiver antenna's orientation vectors. The coeffi- 



cients of the fitting components are preferably functions 
of earth formation parameters such as resistivity of for- 
mation beds, location of the logging instrument, borehole 
deviation, and azimuth angle at the location of the log- 
ging instrument. The fitting coefficients preferably include 
constant, sinO, cosO, sin2cD and cos20 terms that define 
an iterative fitting algorithm useful for determining the 
azimuthal dependence of the directional measurements. 
The measured first and second voltage signals are prefer- 
ably complex voltage signals. Accordingly, in this embod- 
iment, the inventive method further includes the steps of 
calculating the phase-shift and attenuation values from 
the fitting coefficients for the measured first and second 
voltage signals, and combining the calculated phase-shift 
and attenuation values for the measured first and second 
voltage signals to generate a symmetrized or anti- 
symmetrized measurement. The phase-shift and attenua- 
tion values can be obtained by taking the logarithm of the 
ratio of the complex voltage signals obtained from the fit- 
ting expression at two azimuthal angles, preferably az- 
imuthal angles that are 0 and 180 degrees from a deter- 
mined bedding azimuth. 
[0048] | n another aspect, the present invention provides for the 



characterization of noise in the directional measurements. 
Thus, in embodiments wherein only one of the first an- 
tennas has a tilted magnetic dipole, the noise of the mea- 
sured first and second voltage signals can be character- 
ized using the second harmonic coefficients. In embodi- 
ments wherein each of the first antennas has one of a 
tilted and a transverse magnetic dipole, and the fitting co- 
efficients include third harmonic coefficients, the noise of 
the measured first and second voltage signals can be 
characterized using the third harmonic coefficients. In still 
other cases, the noise of the first and second measured 
voltage signals can be characterized by combining the 
first and second measured voltage signals. 
[0049] jhe present invention is also adaptive to conditions 

wherein the logging instrument is not being rotated, such 
as, e.g., during cessation of drill string rotation when di- 
rectionally drilling with a mud motor assembly. Under 
such conditions, the second antennas are preferably sym- 
metric with respect to the first antennas. The azimuth of a 
formation bed of interest is determined by combining the 
first and second antenna couplings (described above), and 
determining the constant and first harmonic coefficients 
from the measured first and second voltage signals ac- 



quired when the instrument is not rotating. The deter- 
mined coefficients can then be used to execute the fitting 
step when the instrument is again rotating. 
[0050] | n another aspect, the present invention provides a 

method and system for using the extracted directional 
signals to obtain boundary distances and make geosteer- 
ing decisions. The extracted directional signals are used 
to extract the distances to the boundary with two tech- 
niques. A simple cross-plot is used for simple one- 
boundary-scenario or known resistivity profile of the for- 
mation, while inversion techniques are used for more 
complex situations and to build consistent structure mod- 
els. A two-dimensional cross plot encompasses measure- 
ments that are sensitive to resistivity and a measurement 
that is sensitive to distance, if shoulder resistivity is 
known. The alternative is to use cross-plot of two direc- 
tional measurements to get the distance and resistivity of 
the shoulder for known bed resistivity. A three-di- 
mensional crossplot can also be generated for one 
boundary situation where each of the formation resistivity, 
shoulder resistivity, and distance to the boundary can be 
generated. The inversion technique uses multiple mea- 
surements to invert for different models and find the best 



match. The inversion can be done downhole or at the sur- 
face. 

[0051] a Graphical User Interface (GUI) is part of the system for 
visualization of the formation model as well as the mea- 
surements and inversion results. The GUI facilitates inter- 
active definition of inversion parameters, selection of 
measurements and models for improving the interpreta- 
tion and generating consistent structure models.The iter- 
ative fitting algorithm preferably includes the steps of: 

initialize P 0 and Z7 0 ; 

for m = \ to Nsamples 

P • r T • r • P 

P y P — JW-1 JW-1 M-1 W-l 

next m; 
return (U); 

[0052] where.Nsamples is the total number of samples acquired in 
one cycle, Mis the dimension of the approximate function 
vector (number of approximation functions), U is the vec- 
tor of fitting coefficients of dimension a matrix of dimen- 
sion M, r is the vector of approximate function values at 
each measure position of dimension M, and P is a matrix 



of dimension Mx M. 

[0053] | n one embodiment, the iterative fitting algorithm deter- 
mines if the fit error is below a pre-defined threshold, and 
if U converges to a value that is representative of the fit- 
ting coefficients. 

[0054] in one embodiment, the iterative fitting algorithm employs 
an integer implementation. Preferably, the integer imple- 
mentation is employed when the logging instrument is 
azimuthally-rotated at relatively high rates, and a signifi- 
cant number of response channels require fitting. 

[0055] in one embodiment, the fitting coefficients are used to 
determine the orientation of a formation bed. The mea- 
sured first and second voltage signals are preferably com- 
plex voltage signals. The orientation of the formation bed 
with respect to the azimuthal angle reference for each 
channel of directional measurement is determined ac- 
cording to: 




= tan l [ 



] 




where C is the real or imaginary part of the coefficient of 



sinO, and C^is the coefficient of cosa from the fitting. A 
common azimuthal angle for the first and second voltage 
signals can be calculated using weighted averaging of the 
fitting coefficients for real and imaginary parts of the 
measured voltage signals. The amplitude and phase of the 
measured voltage signal can be calculated at an assumed 
normal direction to a bed boundary of interest. Phase shift 
and attenuation can be determined by taking propagation 
measurements for two azimuth angles, e.g., O and O 

a a bed bed 

+ 180°. The signals from the fitting coefficients for the 
first and second voltage measurements are preferably 
combined to produce signals necessary to determine the 
distance to bed boundaries of interest. 
[0056] other aspects of the present invention relate to formation 
characterization through the use of crossplots. One 
method includes the step of cross-plotting two directional 
logging measurements acquired from an instrument dis- 
posed in a borehole intersecting the formation to obtain a 
distance to at least one formation boundary and a resis- 
tivity for at least one formation bed. The cross-plotting is 
achieved using a one-boundary model. The obtained re- 
sistivity is the shoulder-bed resistivity, and the obtained 
distance is the closest distance to the shoulder-bed. 



[0057] Another method includes the step of cross-plotting a re- 
sistivity and a directional measurement determined using 
an instrument disposed in a borehole intersecting the for- 
mation to obtain a distance to at least one formation 
boundary and a resistivity for at least one formation bed. 
The cross-plotting is achieved using a one-boundary 
model. The obtained resistivity is the shoulder-bed resis- 
tivity, and the obtained distance is the closest distance to 
the shoulder-bed. 

[0058] a further method includes the step of cross-plotting a re- 
sistivity and two directional measurements determined 
using an instrument disposed in a borehole intersecting 
the formation to obtain a distance to at least one forma- 
tion boundary and a resistivity for at least two formation 
beds. The cross-plotting is achieved using a one- 
boundary model. The obtained resistivities are the bed 
and shoulder-bed resistivities, and the obtained distance 
is the closest distance to the shoulder-bed. The deter- 
mined boundary distance and bed resistivity can be used 
to make drilling decisions. 

[0059] Another aspect of the present invention relates to the use 
of an inversion technique to interpret the directional mea- 
surements for geosteering applications. Here, the iterative 



fitting algorithm is useful for selected real-time direc- 
tional measurements having utility in geosteering. An ap- 
propriate inversion model is selected for the selected 
real-time directional measurements. Once the selected 
model is verified to be consistent with other information, 
it is used to make drilling decisions. 

[0060] The model-selection step preferably includes running a 
plurality of models comprising model-types 
of:homogenous isotropic (single parameter: resistiv- 
ity);homogenous anisotropic (two parameters: Rh and 
Rv);single boundary isotropic formation, boundary above 
or below (three parameters: Rbed, Rshoulder and distance 
to boundary);single boundary anisotropic formation, 
boundary above or below (four parameters: Rbed_h, 
Rbed_v, Rshoulder and distance to boundary);two bound- 
ary isotropic formation three parameters: (five parame- 
ters: Rbed, Rshoulder.up, Rshoulder.down and distance 
to boundary above and below the tool); andtwo boundary 
anisotropic formation three parameters: (six parameters: 
Rbed_h, Rbed_v, Rshoulder_up, Rshoulder.down and dis- 
tance to boundary above and below the tool). 

[0061] The model selection step preferably further includes cre- 
ating a visualization of the selected directional measure- 



merits. 

[0062] | n a particular embodiment, the model-selection step in- 
cludes identifying known formation parameters, interac- 
tively choosing the models with which to invert the se- 
lected directional measurements, and selecting the sim- 
plest model that fits the known information. 

[0063] The model-verifying step includes comparing the selected 
model to known geological characteristics and other mea- 
sured formation parameters, and updating the selected 
model if the selected model is not consistent with the 
known information. 

[0064] The updating step includes refining the selected model 
based upon one of trends, prior knowledge, external in- 
formation, and a combination thereof. Appropriate inver- 
sion parameters are selected, and ranges for the selected 
parameters are defined. The inversion model is preferably 
updated by adding more formation beds. Some of the 
real-time directional measurements can be re-weighted 
or eliminated, and the resulting real-time directional 
measurements re-inverted to the updated model. 

[0065] The crossplotting methods according to the present in- 
vention preferably include defining an appropriate model, 
selecting appropriate directional measurements, inputting 



the selected measurements to the defined model to gen- 
erate the cross-plot, and generating a visual representa- 
tion of the cross-plot. The crossplot can be updated with 
further real-time measurements. 
[0066] a still further aspect of the present invention relates to an 
apparatus for measuring characteristics of earth forma- 
tions surrounding a borehole. The apparatus includes a 
logging instrument adapted for disposal within the bore- 
hole. The logging instrument has a longitudinal axis and 
is equipped with first and second transmitter-receiver an- 
tenna pairs. The first transmitter-receiver antenna pair in- 
cludes a first transmitter antenna having a magnetic 
dipole oriented in a first direction with respect to the lon- 
gitudinal axis of the logging instrument, and a first re- 
ceiver antenna located a first distance away from the first 
transmitter antenna and having a magnetic dipole ori- 
ented in a second direction, the first and second direc- 
tions being different. The magnetic dipoles of the first 
transmitter and receiver antennas define a plane that in- 
cludes the longitudinal axis of the logging instrument. 
The second transmitter-receiver antenna pair include a 
second transmitter antenna having a magnetic dipole ori- 
ented in the second direction with respect to the longitu- 



dinal axis of the logging instrument, and a second re- 
ceiver antenna located the first distance away from the 
second transmitter antenna and having a magnetic dipole 
oriented in the first direction. The magnetic dipoles of the 
second transmitter and receiver antennas define a plane 
that includes the longitudinal axis of the logging instru- 
ment. The apparatus further includes a toolface sensor for 
continuously indicating the azimuthal orientation of the 
logging instrument, and a controller for controlling the 
first and second transmitter-receiver antennas pairs so as 
to selectively transmit electromagnetic energy into the 
formation and measure the voltage signals associated 
with the transmitted electromagnetic energy as a function 
of the azimuthal orientation of the logging instrument. 

[0067] | n a particular embodiment, the second transmitter-re- 
ceiver antenna pairs are oriented at a first azimuthal angle 
(e.g., 90 degrees) with respect to the first transmitter-re- 
ceiver antenna pairs about the longitudinal axis of the 
logging instrument. 

[0068] | n various embodiments, the measured formation charac- 
teristics include resistivity, and earth formation geometry 
information such as dip, azimuth, and bed thickness. 

[0069] The first and second directions are variable, and can, e.g., 



be substantially collinear with the longitudinal axis of the 
logging instrument or substantially 45 degrees from the 
longitudinal axis of the logging instrument. 

[0070] | n a particular embodiment, the first and second transmit- 
ter-receiver antenna pairs are located at the same physi- 
cal positions on the logging instrument. 

[0071] | n a particular embodiment, each of the transmitters and 
receivers have transceiver capabilities. 

[0072] The toolface sensor can employ magnetometers to indi- 
cate the azimuthal orientation of the logging instrument 
with respect to earth"s magnetic north, or gravitation sen- 
sors to indicate the azimuthal orientation of the logging 
instrument with respect to the earth"s gravity vector. 

[0073] jhe inventive apparatus preferably further includes a CPU 
for processing the measured voltage signals within the 
borehole, a telemetry apparatus for transmitting the mea- 
sured signals and CPU-processed results from the bore- 
hole to the surface, and a surface system for further pro- 
cessing measured signals together with other measure- 
ments to generate and display selected parameters of a 
consistent earth model. 
Brief Description of Drawings 



[0074] so that the above recited features and advantages of the 



present invention can be understood in detail, a more 
particular description of the invention, briefly summarized 
above, may be had by reference to the embodiments 
thereof that are illustrated in the appended drawings. It is 
to be noted, however, that the appended drawings illus- 
trate only typical embodiments of this invention and are 
therefore not to be considered limiting of its scope, for 
the invention may admit to other equally effective embod- 
iments. 

[0075] Figure 1 shows schematic diagrams of prior art induction 
or propagation tools. 

[0076] Figures 2A and 2B are elevational views showing eddy 

currents induced by a logging tool in a borehole penetrat- 
ing a formation without and with a relative dip, respec- 
tively. 

[0077] Figure 3 is an elevational representation of a conventional 
rotary drilling string in which the present invention may 
be employed to advantage. 

[0078] Figure 4 is a schematic representation of a basic direc- 
tional measurement logging tool having symmetrical 
transmitter and receiver antenna pairs. 

[0079] Figure 5A is a schematic representation of a directional 

measurement logging tool having a TRR configuration that 



is insensitive to anisotropy at any dip angle, in accordance 
with one aspect of the present invention. 

[0080] Figure 5B shows plots of the directional propagation re- 
sponse for a three-layer formation using a logging tool 
according to Figure 5A. 

[0081] Figure 6 shows a graphical comparison of a binning tech- 
nique and an on-the-fly fitting technique according to 
one aspect of the present invention. 

[0082] Figure 7 shows plots representing the convergence and 
error in a target coefficient set resulting from a fitting 
technique according to one aspect of the present inven- 
tion, as implemented through an integer algorithm in a 
digital signal processor in accordance with another aspect 
of the present invention. 

[0083] Figure 8 shows a cross-plot chart representing a conver- 
sional resistivity measurement and a directional measure- 
ment as used to obtain the resistivity of a formation bed 
and the distance to the bed boundaries, for a known 
shoulder-bed resistivity. 

[0084] Figure 9 shows the results of an inversion technique ap- 
plied to obtain the resistivity and point-by-point bound- 
ary locations of a formation bed. 

[0085] Figure 10 shows a cross-plot-based inversion technique 



applied to interpret the directional resistivity measure- 
ments. 

[0086] Figure 11 is a work flow diagram for geosteering in accor- 
dance with one aspect of the present invention. 

[0087] Figure 12 is a work flow diagram for a cross-plot having 
utility for determining bed boundary distances and apply- 
ing such distances in real-time geosteering, in accordance 
with further aspects of the present invention. 

[0088] Figure 13 shows a block diagram representing the struc- 
ture of a GUI for displaying inversion-based formation 
bed models. 

[0089] Figure 14 shows a computer-generated visualization of 
geosteering inversion associated with the GUI of Figure 
13. 

Detailed Description 

[0090] Figure 3 illustrates a conventional drilling rig and drill 
string in which the present invention can be utilized to 
advantage. A land-based platform and derrick assembly 
10 are positioned over a wellbore 11 penetrating a sub- 
surface formation F. In the illustrated embodiment, the 
wellbore 11 is formed by rotary drilling in a manner that is 
well known. Those of ordinary skill in the art given the 
benefit of this disclosure will appreciate, however, that the 



present invention also finds application in directional 
drilling applications as well as rotary drilling, and is not 
limited to land-based rigs. 

[0091] a drill string 12 is suspended within the wellbore 11 and 
includes a drill bit 15 at its lower end. The drill string 12 
is rotated by a rotary table 16, energized by means not 
shown, which engages a kelly 17 at the upper end of the 
drill string. The drill string 12 is suspended from a hook 
18, attached to a traveling block (also not shown), 
through the kelly 17 and a rotary swivel 19 which permits 
rotation of the drill string relative to the hook. 

[0092] Drilling fluid or mud 26 is stored in a pit 27 formed at the 
well site. A pump 29 delivers the drilling fluid 26 to the 
interior of the drill string 12 via a port in the swivel 19, 
inducing the drilling fluid to flow downwardly through the 
drill string 12 as indicated by the directional arrow 9. The 
drilling fluid exits the drill string 12 via ports in the drill 
bit 15, and then circulates upwardly through the region 
between the outside of the drill string and the wall of the 
wellbore, called the annulus, as indicated by the direction 
arrows 32. In this manner, the drilling fluid lubricates the 
drill bit 15 and carries formation cuttings up to the sur- 
face as it is returned to the pit 27 for recirculation. 



[0093] The drill string 12 further includes a bottomhole assem- 
bly, generally referred to as 34, near the drill bit 15 (in 
other words, within several drill collar lengths from the 
drill bit). The bottomhole assembly includes capabilities 
for measuring, processing, and storing information, as 
well as communicating with the surface. The bottomhole 
assembly 34 thus includes, among other things, a mea- 
suring and local communications apparatus 36 for deter- 
mining and communicating the resistivity of the formation 
F surrounding the wellbore 11. The communications ap- 
paratus 36, also known as a resistivity tool, includes a 
first pair of transmitting / receiving antennas T, R, as well 
as a second pair of transmitting / receiving antennas T", 
R". The second pair of antennas T", R" are symmetric with 
respect to the first pair of antennas T, R, as is described in 
greater detail below. The resistivity tool 36 further in- 
cludes a controller to control the acquisition of data, as is 
known in the art. 

[0094] The BHA 34 further includes instruments housed within 

drill collars 38, 39 for performing various other measure- 
ment functions, such as measurement of the natural radi- 
ation, density (gamma ray or neutron), and pore pressure 
of the formation F. At least some of the drill collars are 



equipped with stabilizers 37, as are well known in the art. 
[0095] a surface/local communications subassembly 40 is also 
included in the BHA 34, just above the drill collar 39. The 
subassembly 40 includes a toroidal antenna 42 used for 
local communication with the resistivity tool 36 (although 
other known local-communication means may be em- 
ployed to advantage), and a known type of acoustic 
telemetry system that communicates with a similar system 
(not shown) at the earth's surface via signals carried in the 
drilling fluid or mud. Thus, the telemetry system in the 
subassembly 40 includes an acoustic transmitter that 
generates an acoustic signal in the drilling fluid (a.k.a., 
"mud-pulse") that is representative of measured downhole 
parameters. 

[0096] The generated acoustical signal is received at the surface 
by transducers represented by reference numeral 31. The 
transducers, for example, piezoelectric transducers, con- 
vert the received acoustical signals to electronic signals. 
The output of the transducers 31 is coupled to an uphole 
receiving subsystem 90, which demodulates the transmit- 
ted signals. The output of the receiving subsystem 90 is 
then coupled to a computer processor 85 and a recorder 
45. The processor 85 may be used to determine the for- 



mation resistivity profile (among other things) on a "real 
time" basis while logging or subsequently by accessing 
the recorded data from the recorder 45. The computer 
processor is coupled to a monitor 92 that employs a 
graphical user interface ("GUI") through which the mea- 
sured downhole parameters and particular results derived 
therefrom (e.g., resistivity profiles) are graphically pre- 
sented to a user. 

[0097] A n uphole transmitting system 95 is also provided fro re- 
ceiving input commands from the user (e.g., via the GUI in 
monitor 92), and is operative to selectively interrupt the 
operation of the pump 29 in a manner that is detectable 
by transducers 99 in the subassembly 40. In this manner, 
there is two-way communication between the subassem- 
bly 40 and the uphole equipment. A suitable subassembly 
40 is described in greater detail in U.S. Patents Nos. 
5,235,285 and 5,517,464, both of which are assigned to 
the assignee of the present invention. Those skilled in the 
art will appreciate that alternative acoustic techniques, as 
well as other telemetry means (e.g., electromechanical, 
electromagnetic), can be employed for communication 
with the surface. 

[0098] Azimuthal Dependence of the Directional Measurement 



and the New Propagation-style Equivalent 

[0099] t wo types of coil antennas are used to compose measure- 
ments with directional sensitivity. One type achieves its 
directional sensitivity by having the antenna either offset, 
e.g., from the center of a logging tool"s longitudinal axis, 
or partially covered. Directional measurements can also be 
made with an antenna configured so that its magnetic 
moment is not aligned with the longitudinal axis of the 
tool carrying the antenna. The present invention relates to 
the second type of directionally-sensitive antenna. 

[0100] Figure 4 schematically illustrates a basic resistivity tool 36 
for directional electromagnetic (EM) wave measurement. 
The tool 36 includes a transmitter antenna T that fires an 
EM wave of some frequency f and a receiver antenna R 
that is some distance L away. Also shown is the symmetric 
pair (T", R") in accordance with the teachings of U.S. 
Patent Application Publication No. 20003/0085707 
("Minerbo et al") assigned to the assignee of the present 
invention. For clarity and simplification, the discussion 
that follows will be limited to the transmitter antenna T 
and the receiver antenna R, although it is generally appli- 
cable to the symmetric antenna pair, T" and R". It should 
be noted that although the tiled moment of the two sym- 



metrization pairs are on the same plane in Figure 4, this is 
not required in the current invention. As will be clear in 
the subsequent discussion, signals from two pairs that 
have their moment in different planes can still be added 
together to achieve equivalent results if the extracted co- 
efficients or directional phase-shift or attenuation are 
used in the symmetrization operation. 
[0101] | n operation, the receiver antenna R will register a voltage 
V induced by the EM wave from the transmitter antenna 

RT 1 

T and its secondary currents produced in the formation 
penetrated by the borehole containing the logging tool 
36. Both antennas T and R are fixed on the tool 36 and 
thus rotate with the tool. The antenna orientations may be 
assumed to form angles 6 T for the transmitter antenna T, 
and 6 for the receiver antenna R. The azimuthal variation 

R 

of the coupling voltage as the tool rotates can then be ex- 
pressed in terms of the coupling of Cartesian components 
of the magnetic dipoles as: 

f„ r (« = W„ cos 6 r cos B t +^)smi3 r «fi3 fl ] 



(1.1) 



where a set of complex coefficients 



c c c c c 



has been defined to represent the amplitudes of the dif- 
ferent components of the measured formation response. 
The complex coefficients are thus defined as: 



[0102] According to one aspect of the present invention, it is rec- 
ognized that these coefficients are functions of formation 
resistivity, borehole deviation, and azimuth angle at the 
tool location. 

[0103] with a symmetrization operation, i.e., 




(1.2) 




, Eq. (1.1) is simplified to: 



— Cj, (0 r , ff,) cosjj* + C X £8 T , 8 X ) sin$ 



(1.3) 

[0104] All the second-order harmonics (C , C 2 ) disappear after 
the subtraction because they are symmetric with respect 
to the exchange of transmitter and receiver tilt angles. 
Thus anti-symmetrization simplifies azimuthal variation 
of the anti-symmetrized signal. 

[0105] At this stage, the reference point of the azimuthal angle is 
arbitrary. For planner geometry, if we choose angle 0 ref- 
erence point as the direction normal to the bedding plane, 
then 




by symmetry and 




would have a pure cosO dependence. In actual application, 
the orientation of the bedding is unknown. However, 



given any reference, the bedding orientation can be calcu- 
lated by: 



. \v -v 

■■tan* A— 

V -V 



(1.4) 

With rotation 

will be normal to the bedding and thus 



fid 



is exactly 



aside from a multiplication constant 



2 sin(fl r - fi £ ) 



[0106] Once the voltage at each of the receiver coils due to each 
of the transmitter coils is determined, the total measure- 



ment can be determined: by adding the voltages in the 
case of an induction tool; or by taking the complex ratio 
of the voltages in the case of a propagation tool. For ex- 
ample, for the propagation logging device of FIG. 4, the 
absolute value of the voltage at each receiver can be ob- 
tained as the square root of the sum of squares of the real 
and imaginary parts of the complex voltage (Eq. 1.1), and 
the ratio of the absolute values provides the attenuation, 
from which the attenuation-determined resistivity R can 
be obtained (resistivity of formations at a relatively deep 
depth of investigation around the receivers). The phase 
for each receiver is obtained from the arc-tangent of the 
ratio of the imaginary and real parts of the complex volt- 
age, and the phase shift is the difference in phase at the 
two receivers. The phase-shift-determined resistivity R 

ps 

can then be obtained (resistivity of formations at a rela- 
tively shallow depth of investigation around the receivers). 
[0107] For propagation-style measurements, the difference of 
the logarithmic of the voltages (or the ratio) between two 
measurements is taken. Following the teachings of 
Minerbo et al, we take the amplitude of the azimuthal re- 
sponse, i.e., the difference in phase-shift and attenuation 
of measurement, at angle <$> and that at (O + 180), evalu- 



ated at the maximum of the voltage response. This leads 
to approximately, from Eqs. (1.1-2 



V R T 1 ^ = C„ (0 r A )^- ',, (-Or ■ 0» ) c " s <l> ^„ (9t , 6 R ) sin 0 4C 2f (8, , 8 K ) cos 14 (8 T , 8 R ) sin 14 
F ftr (IW) + «» "<:,',(()_ .8~~\-C' \o : .O k '\aostt, C u iO r .O, l )smo iC if, |),)«:fi K-,,((> f .0,)Tin2 t » 
= ( C 0 (Or , 0, ) ->C„ (0 r , 0 S ) cos -tC„ (<?,. , ) sin (j, 
~ ' C„ (0 r ,<?„) «? r A)«« 2«i ^(^Jsin 2<* 

sl + 2 , = 1 

r zz cosO J .wsO ll + Jl^+r w ] S mO T sinO R - ( .- I ^ J sin 8 T smQ R sin 2(4 - -WV^ - I'^Jsin^.™^ cos2«> 



(1.5) 



The maximum of 



V 



is achieved at 0= 0 if x is chosen to be the direction nor- 
mal to the bedding. Evaluated at the angle 0 = 0, Eq. (1.5) 
produces: 

[V sintf cose + V cose, ' sinff ] 

^ i . [ K ZX 1 K 



This, however, is still not the pure xz-zx type of re- 
sponses that are desired, i.e., which are insensitive to 
bedding anisotropy and dip angle. 
[0108] The present invention relates to directional measurements 
that are insensitive to anisotropy of the formation at a 
wide range of dip angles and over a wide frequency range. 
As mentioned above, particular embodiments of the in- 
vention are based on anti-symmetrized antenna configu- 
rations or systems. Now with a symmetrization procedure 



ZZ 1 K XX 1 K 




as prescribed by Minerbo et al, we have: 



V^AA) V^AA) ~ ^ C0S V 0S V^ an V^R 

Ci-7) 

This again is similar to the response of the induction type, 
although the denominator still has some components that 
are not simply [xz-zx]. This proves that the symmetriza- 
tion procedure for propagation style measurement can 
produce responses similar to that of the symmetrized in- 
duction type, but not a pure type. It is also true that prop- 
agation measurement can be done at two arbitrary orien- 
tations in the azimuthal response. 
[0109] The final response of this analysis contains components 

from different types of couplings, which at first seems un- 
desirable in light of conventional logging techniques. 
However, this response enables an improved method of 
making the propagation measurement that is simple and 
is much closer to the induction type. The inventive direc- 
tional measurement is achieved by observing that differ- 
ent types of couplings naturally separate into different az- 
imuthal dependences. Note that we can use C , C , C to 

K 0 lc Is 

make a cleaner directional measurement. 
[0110] a determined bedding orientation O may be assumed 

bed 



through Eq. (1.4), such that all of the angles 

^ and x,j? 

are referenced through this direction. In this case, 
(1.1) simplifies to: 

+[V xc sin^ cos^ + ^ cos^.i/n^]cos^ 
+ [i(^-^)sin^ S in^]co S 2 < it 
= C, {9 T , 9 K ) -H^ {S T , B K ) co 5^ (0 r , 0, ) co s 2$ 
(1.3) 

Then a propagation measurement can be simply defined through 

C t (8 rt 8 H ) -C^Sr^J 
(1.9) 

and 

Nowvath symmetriiation, the first egression is obtained as: 
C\(8 r , 8 V ) m C t (8 K ,8^ +Cj£8 x ,8 r ) 

J 7 costf^ cos^ + -[V +V ]dn8 IT ,siyi8 r , 
zz T R 2 xx JF T R 

(1.10) 



which is only slightly different from Eq. (1.7). 



Thus, for both induction and propagation type of twocoil 
(TR) measurements, the analysis of the different compo- 
nent should be done on the voltage V rt (O) level. This 
produces the exact sinO, cos*, sin2<D, and cos20 behav- 
iors that can be extracted through a novel-processing al- 
gorithm described in the next section. 
[° 112 ] It should again be noted that in special cases when either 

8 T - 0 or 8 K - 0 , 

both the sin20 and cos2c|> terms vanish. The voltage de- 
pendence is simply based upon the sinO and cosO terms. 

[0113] one important aspect of the propagation style measure- 
ments of phase-shift and attenuation is that it is naturally 
suited for "while drilling" measurements, for which de- 
tailed characterization of thermal electronics drift under 
downhole conditions is difficult to achieve. The directional 
phase-shift and attenuation measurements defined herein 
have the benefit of a traditional borehole-compensated 
propagation resistivity tool: the transmitter and receiver 
antenna characteristic and the drift of the receiver elec- 
tronics are all dropped out of the measurement. 

[° 114 ] It should also be noted that in this invention the sym- 

metrization process is done with the extracted coefficients 
which can be processed independently of the actual rela- 



tive azimuth of the two TR pairs. If the TR pair orientation 
is physically rotated to an additional angle 0 Q around the 
tool axis, the response is described with exactly the same 
expression except that O is replaced by 

[0115] This has important implications to the actual coil configu- 
rations. By using two pairs of TR measurements, but with 
different azimuthal orientations, the symmetrization pro- 
cess during rotation can be done exactly as if the two 
pairs are on the same plane. However, when the tool is 
sliding, the measured signal from the two orientations can 
be used to construct the required directional measure- 
ments by simply assuming the azimuth of the tool and the 
bedding did not change from the last time when the tool 
is still rotating. Thus, we can use Eq. 1.3 to obtain the 
amplitude Clc and Cls required to construct the propa- 
gation measurements. 

[0116] These analyses can be extended straightforwardly to the 
traditional TRR type of measurements, as described in 
Minerbo et al. One skilled in the art can easily show that 
this procedure produces essentially the same response as 
indicated above, but with twice the signal when the spac- 



ing between the receiver pair is much smaller comparing 
with the TR spacing. The directional signals from the two 
receivers simply add. 
[01 1 7] Figure 5A shows a TRR configuration that is insensitive to 
anisotropy at any dip angle, and Figure 5B shows re- 
sponses according to this configuration. The transmitter 
antenna Tl is energized and the phase shift and attenua- 
tion from the receiver antennas Rll, R12 is measured. 
Then, the transmitter antenna T2 is energized and the 
phase shift and attenuation from the receiver antennas 
R21, R22 is measured. The tool reading corresponds to 
the differences between these two sets of measurements. 
Since the individual measurements are identical in a ho- 
mogeneous medium at any angle and with any anisotropy, 
the tool readings is zero in a homogenous medium at any 
dip. 

[0118] The measurement responses in a three-layer anisotropic 
formation are shown in Figure 5B. The tool reading is zero 
far from the boundary at any dip, and there is little sensi- 
tivity to anisotropy close to the boundary. Separation in 
responses comes from the fact that propagation re- 
sponses are not symmetric if the transmitter and receiver 
location are interchanged. Making an up/down measure- 



ment contains only the directional information, even close 
to the boundary. It should be observed that attenuation 
responses are practically overlapping for different dip if all 
antennas are in the same medium, similarly to ideal XZ-ZX 
induction measurement (described above). The phase shift 
measurements are also overlapping, although responses 
are double-valued in the conductive bed (1 S/m). 

[° 119 ] Digital Signal Processor (DSP) Algorithms to Extract the 
Signal from the Azimuthal Dependence 

[0120] it j S c | ear f rom the previous analysis that the best way to 
extract the directional measurements is to analyze the 
voltage signals and extract the relevant measurements 
and then synthesizes them with symmetrization for final 
use with geosteering. The traditional method to deal with 
azimuthal data is to bin them into small bins of equal az- 
imuthal span (see left portion of Figure 6), as was done 
with resistivity or density borehole images. Then, the up 
and down bins could be identified by locating the maxi- 
mum and minimum. However, there are many disadvan- 
tages of this technique: 

[0121] i. binning reduces the angle resolution of the measure- 
ment; 

[0122] 2. binning results in wasted data since none of the data 



outside the up and down bins is used for distance- 
to-boundary calculation; 

[° 123 ] 3. in stick and slip occurrences, binning is not uniform and 
thus the up and down bins may be empty or have small 
samples, thus inducing error in the up/down measure- 
ments and possibly producing the wrong identification of 
peak values in the worst case; 

[0124] 4. binning consumes more memory. 

[0125] one aspect of the present invention is principally con- 
cerned with determining the coefficients of the sinO, cosO, 
sin20 and cos20 terms, which define a new algorithm to 
be used instead of conventional binning. This inventive 
technique is referred to as "on-the-fly"fitting of the az- 
imuthal response according to extracts of the relevant sin 
and cos terms of the directional measurements, taken it- 
eratively (see plotted points of Figure 9). Such a fitting al- 
gorithm is done in a DSP through an integer algorithm so 
it is fast enough to be performed for all channels within 
the 4-ms of sampling time. The precise use of azimuth 
angle information, and the randomization of the acquisi- 
tion sequences, makes the algorithm robust to tolerate ir- 
regular tool rotation as well as stick-and-slip under rough 
drilling conditions. This way, all the data are used to ob- 



tain the up/down signal instead of only the data in the 
two bins, thus improving the signal-to-noise ratio in the 
measurement. The use of precise azimuth angles also 
makes the determined bedding orientation more precise. 
[0126] The detail algorithm can be described as follows. 

[0127] Floating point implementation: starting with an initial 

value of matrix P q and vector U o , then proceeding to the 
algorithm described below (also represented graphically in 
the right-hand portion of Figure 6) with measurement 

y($i) and basis r = (j cos£ sin£ cos2£ sin2^) r , where P is a matrix of dimension 
M x M and U andr are vectors of dimension M . M is the dimension of the basis 
function. After iteration!/ , thenU 

will converge to a value which represents the coefficients 
of the expression. This algorithm is stable and conver- 
gence is usually achieved within 10-15 iterations.The de- 
tailed algorithm is shown below: 



initialize P () and U f) ; 

for m = 1 to NsampSes 

next m- 
?eturn(Uy, 

vdiere: 

Samples is the total number □ f samples acquired in one cycle, 

is the dimension of the approximate function vector (number of approximation 
functions), ^ 

Uisthevector of fitting coefficients of dimension , M 
7 is the vector of approximate function values at each measure position of dimension 
and M x M 

P is a matrix of dimension 

[0128] integer Algorithms in DSPIn many cases, floating point 
implementation will be too expensive to perform with 
presently available downhole cpu"s because there may be 
hundreds of channels to be fitted and the data acquisition 
for each azimuth angle has to be quite short (ms) in order 
for the angle to be accurate at higher rotation speed. In 
this situation, an integer implementation can be applied, 
with some modification, to improve accuracy (e.g., use 
32-bits for multiplication), perform rescaling to avoid 
overflow, and to accelerate convergence. The values of the 
basis function can also be pre-generated and stored in 
memory so as to be interpolated later to obtain the value 
for the true angle O . 



[° 129 ] In summary, the integer algorithms are: 

[0130] -based on 16-bit integer data representation with 32-bit 

integer intermediate variables; 
[0131] -more efficient but less precise; 

[0132] -adaptive to strategies for improving accuracy and speed 
of convergence: scaling, initialization and re-initialization; 

[0133] -adaptive to strategies to adapt to DSP environment: divi- 
sion, rounding. 

[0134] a convergence flow of the integer version of the algorithm 
is shown in Figure 7. It should be noted that the fitting is 
quite accurate with errors generally less than 1%. A very 
important point of the fitting is that almost all the data 
are used to derive the coefficients (directional measure- 
ment signals), thus improving significantly the signal- 
to-noise ratio. For example, if 32 bins are used, then only 
1/16 of data are used if the up/down binning is imple- 
mented. Using on-the-fly fitting, in contrast, almost all 
the data are taking into account (aside from the initial 
convergence part). 

[0135] since only the relevant signals are extracted in the fitting 
technique (described above), only the useful coefficients 
need to be saved. Thus, in this case, it"s only necessary to 



save 5 coefficients, as compared to 32 if one were to bin 
all the data using the 32-bin example. Those skilled in 
the art will appreciate the advantages of the inventive 
technique which include the accuracy of the extracted sig- 
nal and a particular improvement in the accuracy of the 
azimuthal angle. 

[0136] | n another aspect, the present invention provides for the 
characterization of noise in the directional measurements. 
Thus, in embodiments wherein only one of the first an- 
tenna pair (TR) has a tilted magnetic dipole, the noise of 
the measured voltage signals can be characterized using 
the second harmonic coefficients. In embodiments 
wherein each of the first antennas has either a tilted or a 
transverse magnetic dipole, and the fitting coefficients in- 
clude third harmonic coefficients, the noise of the mea- 
sured voltage signals can be characterized using the third 
harmonic coefficients. In still other cases, the noise of the 
measured voltage signals can be characterized by com- 
bining the signals. 

[0137] The present invention is also adaptive to conditions 

wherein the logging instrument is not being rotated, such 
as, e.g., during cessation of drill string rotation when di- 
rectionally drilling with a mud motor assembly. Under 



such conditions, the second antenna pair (T'R") is prefer- 
ably symmetric with respect to the first antenna pair. The 
azimuth of a formation bed of interest is determined by 
combining these antenna couplings (as described above), 
and determining the constant and first harmonic coeffi- 
cients from the measured voltage signals acquired when 
the instrument is not rotating. The determined coeffi- 
cients can then be used to execute the iterative fitting 
when the logging instrument is again rotating. 
[0138] Boundary Distance for Geosteering UseAccording to an- 
other aspect of the present invention, the determined co- 
efficients can be used to obtain the orientation of the for- 
mation bedding. For each channel of directional measure- 
ments, the orientation of the bedding can be determined 
through Eq. 1.4. It may be further noted that the average 
value from many channels, weighted by a function of the 
relative signal strength of each channel that comprise a 
measurement, can also be used to increase accuracy since 
this orientation should be the same for all channels. The 
amplitude and the phase of the voltage signal 




at the assumed normal direction to the bed can then be 
calculated. 

[0139] The directional phase-shift and attenuation can be calcu- 
lated with Eq. 1.10. Symmetrization is then performed to 
produce the final signal needed to find the distance to the 
boundary. It should be appreciated, however, that the or- 
der of these steps can be varied with similar or identical 
responses. 

[0140] jo obtain the boundary distance in real-time, two tech- 
niques can be employed. For simple models (one bound- 
ary only), a cross-plot of two directional measurements 
enables us to get both the distance to the boundary and 
the formation resistivity of one of the beds. A representa- 
tive cross-plot chart is shown in Figure 8, using the re- 
sponse of an 84 inch, 100kHz attenuation versus a 28 
inch, 2MHz phase shift resistivity (shoulder resistivity be- 
ing R = 0.8Qm, R = 3.Qm). Here, symmetrization effec- 

h v 

tively removes other parameters such as anisotropy and 
dipping from the consideration. The charts for sym- 
metrized measurement are simple. The use of different 
combinations of pair measurements to obtain a consistent 
picture will increase the confidence of the interpretation. 
[0141] Figure 9 illustrates the use of cross-plot based inversion 



to interpret the directional measurements. Two directional 
84 inch 400 kHz measurements are used to interpret the 
measurements using a single boundary model with fixed 
bed resistivity at 1000m, and varied shoulder bed resis- 
tivity (below the tool) and distance. The response is not 
sensitive to bed resistivity, and signal is primarily deter- 
mined by the shoulder-bed resistivity and distance. The 
dots plotted on the captured screen image indicate shoul- 
der-bed resistivity measurements and distance readings. 
The value of resistivity and distance are read from the 
cross-plot and output on the screen. 
[0142] For the one-bed boundary case where the resistivities of 
the beds on both sides of the boundary are known, one 
can use three inputs, at least one of which is a measure- 
ment of the resistivity surrounding the tool (e.g., tool 36 
in Figure 3). 

[0143] For more complex models involving more boundaries, an 
inversion program based on the disclosed technique of 
U.S. Patent No. 6,594,584 is employed. A trajectory of the 
borehole is projected onto an initial model of the forma- 
tion. A logging tool is disposed within a segment of the 
borehole trajectory, and its responses along the segment 
are measured. The expected responses of the tool, ac- 



cording to the model, are also determined. Differences re- 
sulting from a comparison between the expected and the 
measured responses along the segment are then used to 
adjust the model, and the cycle of comparing and adjust- 
ing are repeated until the differences fall below a selected 
threshold. The robustness of the inversion model is im- 
proved with multiple starting points and physical criteria 
to distinguish between solutions. 

[0144] A n inversion algorithm is then applied. This algorithm ac- 
cepts any measurements as input, and then finds the 
most consistent model to the data. Multiple models are 
run, and the best model is selected automatically. A typi- 
cal resistivity profile reconstruction is shown in Figure 
10A, while the corresponding boundary (structure) recon- 
struction is depicted in Figure 10B. These results indicate 
that the solutions are more accurate when approaching 
conductive beds than when approaching resistive beds. 
This is expected since the directional responses are pri- 
marily defined by the conductive bed, and have far better 
sensitivity to resistivity variation in conductive layers 
above and below, then in the resistivity variation of resis- 
tive shoulder beds. 

[0145] Another aspect of the present invention is the utilization 



of the determined boundary distances to make geosteer- 
ing decisions. The overall geosteering workflow can be 
described with reference to Figure 11. First, real-time di- 
rectional measurements are selected at box 110, using 
the on-the-fly fitting technique as described above. The 
directional measurements are then visualized, at box 120 
by showing the raw log responses, or using cross-plot 
charts as shown in Figures 8 and 10 (and described else- 
where herein). 

[0146] Decision point 130 then directs the work flow, depending 
on whether an appropriate model has been identified. If 
the formation layout or some other parameters are 
known, this approach allows the user to fix the known pa- 
rameters. For example the upper shoulder bed resistivity 
or layering may be known to be stable, and it may be 
known with high certainty that there is no boundary be- 
low. Such information allows, e.g., selection of model with 
a single boundary, and fixed shoulder bed resistivity. 
Mathematically, it means that only three parameters need 
to be inverted, horizontal and vertical resistivity of the 
bed (R and R ), and distance to upper boundary. This is 

h v 

an example of full model inversion (box 140), which in- 
sures consistent interpretation, and avoids confusion in 



limiting cases where there may be some non-physical 
model with better mathematical fit of measurement to be 
selected. 

[0147] Alternatively, fast approximated inversions with multiple 
models are run (box 150) along with an automatic model 
selection algorithm (box 160). The algorithm runs differ- 
ent models, from simple (no boundary, isotropic forma- 
tion) to most complex (two distances and anisotropic for- 
mation), including: 

[0148] -homogenous isotropic (single parameter: resistivity); 

[0149] -homogenous anisotropic (two parameters: Rh and Rv); 

[0150] -single boundary isotropic formation, boundary above or 
below (three parameters: Rbed, Rshoulder and distance to 
boundary); 

[0151] -single boundary anisotropic formation, boundary above 
or below (four parameters: Rbed_h, Rbed_v, Rshoulder 
and distance to boundary); 

[0152] . two boundary isotropic formation three parameters: (five 
parameters: Rbed, Rshoulder_up, Rshoulder.down and 
distance to boundary above and below the tool); and 

[0153] -two boundary anisotropic formation three parameters: 
(six parameters: Rbed_h, Rbed_v, Rshoulder_up, Rshoul- 



der_down and distance to boundary above and below the 
tool). 

[0154] a solution is obtained for each model. 

[0155] The model selection algorithm uses physics-based con- 
straints for directional measurements, and imposes con- 
ditions that select the "simplest model that fits the 
data"(Bayes" razor). The classic Akaike Information Crite- 
rion (AIC) or Bayesian information criterion can be used, 
to penalize the model complexity. 

[0156] a model validation is then applied, at decision point 180, 
to determine whether the selected model is consistent 
with prior knowledge about the geological formation or 
other measurements such as gamma-ray, or other mea- 
surements available in the real-time as represented by 
box 170. If the model is consistent with the other data, 
it"s accepted as input (box 190) to the drilling decision 
(box 200), and an appropriate report is generated at box 
210. The implementation of the drilling decision will re- 
sult in further real-time directional measurements (box 
100), which are again input to the fitting and visualization 
steps of boxes 110 and 120. If the selected inversion 
model is inconsistent with the other measurements from 
box 170, then an update of the model is appropriate. In 



this instance, an interactive mode sub-flow 220 is applied 
that is consistent with the complex inversion program of 
the "584 patent. 

[0157] often, the automatic inversion may create interpretation 
that is not consistent. It can be caused by the measure- 
ment noise caused by electronics as well as the "model 
noise", i.e., the fact that the real model is different from 
all the models run in the inversion. For example, borehole, 
invasion, thin layers, cross-bedding and finite tool size 
are not included in models, and they can cause mismatch 
in the fitting process. 

[0158] The flexibility and ability of a user to interactively select 
the common model for a segment of data is key for suc- 
cessful interpretation of the measurements. The software 
defining the interactive mode sub-flow (box 220) has fea- 
tures that allow: model refinement (box 221) based on 
trends, prior knowledge, or an external source of infor- 
mation; constraining or fixing some inversion parameters 
(box 222); removing some measurements that may be 
more affected by the environment not included in the 
model; and re-processing the data (box 223). 

[0159] Figure 12 illustrates a general workflow for creating 
cross-plots of directional measurements that provide 



boundary distances, as described elsewhere herein. The 
creation or modification of a cross-plot is represented by 
sub-flow 20, which is initiated by defining an appropriate 
model (box 30). Appropriate directional measurements 
are then selected according to the ability to determine/pre- 
dict shoulder resistivity (box 40). If shoulder-bed resistiv- 
ity is reliably known as is often the case in big fields 
where many wells are drilled, interpretation is based on 
one resistivity and one directional measurement (box 50), 
to determine the true bed resistivity (shoulder-bed cor- 
rected) and the distance-to-boundary. If shoulder bed re- 
sistivity is not known, it is recommended to use the 
crossplot of two directional measurements (box 60), as il- 
lustrated in Figure 9, to determine the shoulder-bed re- 
sistivity and distance to boundary. The selected measure- 
ments are input to the defined model to generate re- 
sponses, at box 70, which may be visualized (box 80) as 
indicated by Figure 8 and 9. In cases when neither resis- 
tivity is known, but the ranges of their difference is 
known, it is possible to create three-dimensional cross- 
plots, by combining one resistivity and two directional 
measurements. 
[0160] Once the cross-plot has been created or updated 



(sub-flow 20), it may be continuously updated with addi- 
tional real-time measurements (box 100) and visualized 
(box 120"). Processing the cross-plot of directional mea- 
surements yields the distance to one or more formation 
boundaries and the resistivity of one or more beds (box 
220), which may be output and visualized (box 230) to in- 
struct the drilling decision, at box 200. 
[0161] According to another aspect of the present invention, a 

graphical user interface (GUI) has been created to facilitate 
the use of the program and to display the inverted bed- 
ding models for user visualization. A block diagram illus- 
trating the structure of the GUI and different ways of han- 
dling data is shown in Figure 13. Accordingly, the GUI al- 
lows: 

[0162] _ running forward model for given layering 

[0163] _ cross-plots for measurement sensitivity analysis 

[0164] _ inversion of real-time data characterized by: cross- 
plot-based inversion of one resistivity and one distance 
(for given shoulder or bed resistivity); utomatic interpreta- 
tion, using fast inversion based on superposition of single 
boundary responses and automatic model selection, 
physics-based constraints, combined with Akaike Infor- 



mation Criterion or model mixing; interactive interpreta- 
tion based on full inversion of any parameter (up to 6) 
from three-layer model (bed and two shoulders two dis- 
tances and four resistivities), running full layered medium 
model in the inversion loop, allowing model refinements, 
with options to select/deselect some inversion parameters 
and constrain them and/or select and re-weigh available 
measurements; switching between fast and interactive 
mode for model refinement 
[0165] _ visualization of measurement and inversion. 

[0166] An example of an interpretation screen is shown in Figure 
14. Both the measurements that are used in the inversion, 
and the borehole trajectory, are displayed on the screen. 
The inverted results, including the distance to upper and 
lower boundaries, resistivities of bed and two shoulder- 
beds are graphically displayed. The results are shown 
point-by-point, and based on reconstructed positions of 
boundaries. It is possible to infer the stratigraphic dip. It 
should be noted that although the forward model run in 
the inversion loop is a ID-layered medium, the approach 
allows for the building of more complex formation im- 
ages, including 

[0167] -non-parallel boundaries 



[0168] -detection of faults, including sub-seismic faults. 

[0 1 69] Anisotropy Evaluation 

[0170] with A TRR configuration having the transmitter antenna 
and at least one of the receiver antennas being tilted, 
anisotropy measurements can be performed, according to 
the teachings of U.S. Patent Application Publication No. 
2003/020029 ("Omeragic et al"). The detected EM signals 
relating to a phase difference or a magnitude ratio are 
combined to determine the anisotropy. Note that a novel 
borehole compensation technique is disclosed by Omer- 
agic et al for such a measurement. 

[0171] The directional part of such combined measurements can 
be analyzed in the same way as described above. Accord- 
ingly, the component coefficients can be used to derive 
anisotropy through inversion with the technique disclosed 
in U.S. Patent No. 6,594,584 (also described above). This 
is true for all the dip angles. 

[0172] Those skilled in the art will appreciate that in highly devi- 
ated boreholes the traditional propagation measurement 
with axial coils already provides acceptable anisotropy 
sensitivity, obviating the need for the borehole compen- 
sation technique in those boreholes. 



[0173] it will be apparent to those skilled in the art that this in- 
vention may be implemented using one or more suitable 
general-purpose computers having appropriate hardware 
and programmed to perform the processes of the inven- 
tion. The programming may be accomplished through the 
use of one or more program storage devices readable by 
the computer processor and encoding one or more pro- 
grams of instructions executable by the computer for per- 
forming the operations described above. The program 
storage device may take the form of, e.g., one or more 
floppy disks; a CD ROM or other optical disk; a magnetic 
tape; a read-only memory chip (ROM); and other forms of 
the kind well known in the art or subsequently developed. 
The program of instructions may be "object code," i.e., in 
binary form that is executable more-or-less directly by 
the computer; in "source code" that requires compilation 
or interpretation before execution; or in some intermedi- 
ate form such as partially compiled code. The precise 
forms of the program storage device and of the encoding 
of instructions are immaterial here. Thus these processing 
means may be implemented in the surface equipment, in 
the tool, or shared by the two as known in the art. It will 
also be appreciated that the techniques of the invention 



may be used with any type of well logging system, e.g. 
wireline tools, LWD/MWD tools, or LWT tools. 

[0174] it w j|| De understood from the foregoing description that 
various modifications and changes may be made in the 
preferred and alternative embodiments of the present in- 
vention without departing from its true spirit. 

[0175] This description is intended for purposes of illustration 

only and should not be construed in a limiting sense. The 
scope of this invention should be determined only by the 
language of the claims that follow. The term "compris- 
ing"within the claims is intended to mean "including at 
least"such that the recited listing of elements in a claim 
are an open group. "A,""an"and other singular terms are 
intended to include the plural forms thereof unless specif- 
ically excluded. 



